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ABSTRACT 


Through the use of asymptotic and perturbation methods, 
this thesis presents a theoretical study of the flow of a 
Stratified fluid over variable topography as a model of the 
resonant interaction of a submarine's wake with a stratified 
iC tae Be Such resonant interactions may be able to produce 
Significant upstream disturbances. The long time solution 
obtained in our model exhibits growth in time for the resonant 
case, indicating that perhaps nonlinear effects, balanced by 
dispersion, could cause significant upstream disturbances. 
These disturbances could allow the possibility of detecting a 
Submarine by the internal waves its wake generates in the 
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I. INTRODUCTION 


The nonacoustic detection of submarines is an area of 
intensive and well funded research. Although, it is not clear 
whether devices based on nonacoustic means of detection will 
become as useful as those based on acoustic means, the 
military importance of such devices demands pursuing such 
avenues of technological research. One possible nonacoustic 
means of detection is to attempt to detect the submarine 
optically through the use of a laser by examining the way the 
submarine reflects or absorbs blue-green light. Another more 
promising method is to try to detect the submarine by the heat 
it emits. However, an even more promising technique is the 
detection of submarines by the internal waves they generate. 

In the ocean a stable density stratification generally 
exists. As the submarine moves through the density-stratified 
ocean, it leaves behind a wake capable of generating internal 
waves. These internal waves result from the stabilizing 
effect of the stratification: as a water particle is 
displaced, it experiences a buoyant restoring force in the 
opposite direction from which it was displaced. 

Many theoretical and experimental studies have been 
conducted on the subject of the flow of a stratified fluid 
over a localized topography. For the two-dimensional flow of 


an inviscid, incompressible stratified fluid, most studies 


have been confined to steady flow or linearized theory. 
Earlier studies concerned themselves only with a description 
of the downstream stationary lee-wave field. A question that 
has generated some controversy is the phenomenon of upstream 
influence. To what extent do disturbances generated by the 
topography through transient and nonlinear processes propagate 
upstream? 

McIntyre (1972) showed that the only significant upstream 
disturbances were weak second order long-wave motions 
generated by nonlinear interactions in the lee-wave. McIntyre 
(1972) extended the earlier theories of Benjamin (1970) and 
Keady (1971) by performing an expansion using as the small 
parameter the ratio of the height of the topography to the 
total fluid depth. Baines (1977,1979), on the other hand, in 
experimentally examining the flow of a continuously stratified 
fluid over topography, found upstream disturbances that were 
first order in the topography height. These upstream 
disturbances were not only generated by nonlinear processes 
over the topography but also travelled at the long-wave speed. 
Baines also observed that these upstream disturbances were 
very strong when the fluid approached resonance. Here, 
resonance is defined as the near coincidence of the basic flow 
speed and one of the free long-wave speeds. Grimshaw and 
Smyth (1986) presented a theoretical study of stratified fluid 
flow over localized topography for the two-dimensional case 


when the basic flow is near resonance. One interesting 


feature of their results is the appearance of significant 
upstream disturbances as observed by Baines (1977,1979). 
Grimshaw and Smyth showed that solitary waves were formed when 
the localized forcing over the topography produced a positive 
disturbance upstream, and an oscillatory wavetrain was formed 
when the disturbance had negative polarity. 

Motivated by these results, this thesis presents a 
theoretical study of the flow of a stratified fluid over 


localized topography for the three dimensional case. 


II. PROBLEM FORMULATION 


A. DIMENSIONAL PROBLEM 

We consider the three-dimensional flow of an inviscid, 
incompressible stratified fluid. Figure 1 (see page 5) 
provides a sketch of the coordinate system. We assume that 
the basic state has a constant horizontal velocity, V, from 
left to right with equilibrium pressure and density 
respectively as Pp, (z) and tae Here and throughout this 
paper barred variables represent dimensional variables. We 
introduce { as the vertical particle displacement from some 


initial position at time t = 0 to the position at time t =e 


u;+ V+ w= 0, (2 ie) 

p,(2-f) eu + Pe = 0, (2082) 

p, (2-5) BY +p =0, ( 2g) 
p,(2-) DY + p,(Z-$)g + Be = 0, (2.4) 


Figure 1. The Coordinate System 


where 
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Here, u, v, and w are the velocity components relative to the 
basic state which has velocity (V,0,0) with respect to a 
right-handed Cartesian coordinate system. 

The bottom topography is given by z= icy oe where 
initially £ Cy) — 0m For large times the topography 


attains a steady value g(x,y) and hence we assume that 


f(x,y,t) ~ g(x,y) as t7o. (2a 


We also assume that at all times the topography is localized 
where f(x,y,t) > 0 as x+y’ > ». The kinematic bottom boundary 


condition is then 


(258) 
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We are assuming that this "topography" is created by the 
submarine's wake, which when viewed from the moving submarine, 
actually has a steady shape. In this thesis, we investigate 
the interaction of the wake with internal waves. Kallen 
(1987) has investigated this type of interaction ina linear, 
steady state model and comments in detail on the validity of 
assuming that the wake acts as rigid topography. Furthermore, 
from experimental results, Lin and Pao (1979) show the mean 
flow occurs at the wake boundaries implying that part of the 
flow travels around rather than through the wake. It is this 
flow around the wake that generates the internal waves. Our 
model, in fact, does not depend critically on the assumed form 
of "topography" the submarine creates and our purpose here is 
to investigate a time dependent model. 

The upper boundary, assumed to be a free surface, is 
Tocated at z = H + eA) , where H is the vertical distance 
from the surface to the submarine and 1 (xX,y,t) is the free 
surface displacement. One must remember to include not only 
the kinematic boundary condition at the upper boundary but the 
dynamic boundary condition as well. The upper boundary 


conditions take the form: 


t=7 alte z=H+ 4, (Zeca) 


p,(H+n) + p,(H+n) =p,(H) at Z=H+ 7. (2. 9b) 


Here, P, is the perturbation to the pressure D, of the baste 


Staver 


B. DIMENSIONAL SCALING 

The solution of practical problems in fluid mechanics 
requires both theoretical developments and experimental 
results. By grouping significant dimensional quantities into 
dimensionless parameters it is possible to present results in 
a compact form which is applicable to all similar situations. 

In the problem, the characteristic horizontal length scale 
L is determined by the ratio U/N, where U is a characteristic 
horizontal velocity and N is the Brunt-Vaisala frequency at 
the location of the submarine. We shall take the 
characteristic horizontal velocity U to be given by N,H, (U = 
N,H), where N, is the maximum Brunt-Vaisala frequency and H, 
the characteristic vertical length, is the depth of the 


' and H is 200 metersuaae 


submarine. Typically N, is 10° sec 
so U is 2 m/sec. This is a typical order of magnitude 
patrolling speed for a submarine at that depth. Of course, 
O(1) variations in N, and H can be made without changing our 
analysis. 

From the mass conservation equation and with the imposed 
horizontal length and velocity scales, we see that the 
vertical velocity scale, W, is UH/L. The horizontal length 


and velocity scales are also used in determining the 


Characteristic time scale, T, where T is defined by T = L/U = 


1/N. Lastly, the characteristic pressure scale is defined as 
Pp = 9,gH where p, 1S a characteristic density. These scales 
define the parameter 6 = N,H/g, which is a dimensionless 
Measure of the ratio of the buoyancy force to the 


gravitational force. 


C. NONDIMENSIONALIZATION 

By using the various dimensional scales introduced in 
section B, we shall nondimensionalize the equations of motion 
Wee-2.6) along with the boundary conditions (2.8-2.9b). 
However, before deriving the nondimensional equations of 
motion and boundary conditions, we will consider several 
related dimensional equations that will assist in simplifying 
the nondimensional derivation of equations (2.1-2.6) and 
poumdary conditions (2.8-2.9b). We first introduce the 


following nondimensional variables: 
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To nondimensionalize the hydrostatic equation, we begin with 











= -gp,. (2, sige) 
zs g 
Using (2.10) yields 
Dp, (Zz) = Som (2s 
The dimensional Brunt-Vaisala frequency is defined as 
= do 
N) 2 ee (28 188) 
p, az 
where dp,/az < O. Using (2.10), we obtain 
N-H d 
e! ay? Se (2.14) 
g p, az 


ue 6 


N,/H/g, therefore (2.14) simplifies to 





(2. ay 


1G 


In nondimensionalizing the Lagrangian representation of the 


particle's vertical displacement, we start with 


Z=2- €(x,y,z,t). (2.16) 


Meng (2.10) yields 


2 eZ Veneer) (22°07) 


The last equation needed before considering the equations of 


motion and boundary conditions is the dimensional pressure 


equation: 
p(x, y, z, t) ==) GHD). (2) po nua xv, 2, t) C2redeey) 
where P, = p,gHbq(x,y,z,t) and 6 is defined after (2.14). 


Thus, q(x,y,z,t) is the dimensionless perturbation to the 
hydrostatic pressure p,. 

Before deriving the nondimensional equations of motion and 
their associated boundary conditions, we first introduce the 


following nondimensional variables: 


Hak 


ee a ae ee, eee a, 
L L H U 
v= , ee wer , (.. =e 
U U WwW p, 
= = = - (2.19) 
t i. n f 
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where L, H, U, W, p,, and T are the characteristic dimensional 


is a characteristeae 


scales defined in section B and £ h 


amplitude of the topography. 


We begin with the mass conservation equation (2.1). Using 


(2.19) and remembering that W = UH/L, (2.1) becomes 


Ul Vit Wee Oe (Zee) 


To assist in deriving the nondimensional form of the x- 


component of the momentum equations, given by (2.2), we shall 


now derive the nondimensional form for the conveceive 


derivative D/Dt: 


DD. =< a (v+aee st ae ot FS (2 72a) 
Dt Om Ox Oy 0Z 
Using (2.19), (2.21) becomes 


aeZ 


2 =i 7 : 2 (Vou) ' oy : et (22) 
but T = L/U, therefore (2.22) simplifies to 
- = aoe (725) 
where 
ee (Vu) atv twa. (252 4)) 


ieee uting (2.17), (2.18), and (2.23) into (2.2) and using 


(2.19), we obtain 


U Du p,9HB 
ZC) eee a ey ONE 


(2259 


A similar analysis for the y-component of the momentum 


equations (2.3) yields 


ale 


Dv 
p.(Z-$) He +, == 0. (2.26) 


To derive the nondimensional equation for the z-component of 


the momentum equations (2.4), we begin with 


4) Se = 
(2a 0) ==> pe 2-6) Gp eer (2°28) 
p ee $)9 + DP; 

Substituting (2.17), (2.18), and (2.23) into (2.27) angie 


(2.19), we obtain 


W Dw 


pi \zZH ae 7p, 2—().9. + Gp amo. (2.216) 


Using T = L/N,H and W = N,H’/L, (2.28) simplifies to 


2 
WH =) o,(2-0) BY + pj (2-9 + Ie + gBq, = 0. (2 
fe Dt 
Using (2.12) and the definition of 6, B = N/H/g, (2.29) 


Simplifies to 
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Dw 
ep, (2-£) =e + G(o,(2-£) - p(z)) + a = 0, (2.30) 
where 


CZee snl) 


™® 
Ht 
eal ies 


Lastly, we derive the nondimensional vertical velocity 


Meo a1tS definition (2.5). Using (2.19) yields 
ee (32) 
Dike 


(2533) 


om | 
lI 
bh | 
%y 
ct 
N | 
{I 
Hh | 


Using (2.19), we obtain 


¢=af at Zin OE, (2.34) 


where 


ES 


4 2. (2 iy) 
For the free surface kinematic boundary condition, 
esa at z=H+%7 (2.36) 
We use (2.19) to obtain 
¢=yY at Z=1+y%. (2.308 


Next, we nondimensionalize the free surface dynamic boundary 


GOnGIitiIOn in (229n)= 


p,(H+n) + p,(H+n) =p,(H) at 2 = oe (2.38) 
Using (2.19) and (2.18), (2.38) becomes 
Pls yee: = peas at ZzZ=1ty. (2.39) 


To recapitulate our nondimensional versions of equations 
(2.1-2.6) and the surface boundary conditions (2.8-2.9b), the 


nondimensional equations of motion are 
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eee Ve 
fe, Du 
p,(2 ae 
Dv 
Pol Z-S) Te 
2 Dw ali 
ap, (Z- are a Bleo( 2-5) peel 
where € = H/L and 
1B) 0 0 
— = V 
eta tele TE: 
with boundary conditions 
> =arf at 
ae) at 
fet n) + PQ = p (1) ae 


where q@ = h/H. 


+w, = 0, (2.40) 
Ten = 307 (ae) 
pico, (2342) 
+q = 0, (2.43) 
_ Dg 
WwW Fae (2.44) 
+ on + (245) 
Zar (2.46) 
Zz=1+%, (2.47a) 
Zool th, (2 47D) 


An exact analytical solution of these nonlinear partial 


differential equations is not 


available. However, we shall 


ae 


obtain approximate solutions for the case when ¢ defined by 
(2.31) and a@ defined by (2.35) are both small positive 
constants. Since ¢ = H/L = N/N,, the condition ¢ <<3siae 
satisfied when the density changes rapidly across a thin 
thermocline. However, there is some ambiguity in the 
definition of the thermocline. We will consider the 
thermocline to be the region where density variations are 
determined chiefly by temperature variations. This means that 
the submarine is usually submerged to a depth below the 
thermocline. Is this a valid assumption? To avoid acoustic 
detection the submarine submerges to depths below the mixed 
layer or layer depth. Depending on the thickness of the 
thermocline the submarine can be either in the thermocline or 
below it. We shall consider only the case for which the 
submarine is below the thermocline in which the parameter a, 
the ratio of the characteristic amplitude of the topography to 
the submarine's depth defined by (2.35), is also small in this 


Situation. 
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III. THE ASYMPTOTIC EXPANSION OF THE DEPENDENT VARIABLES 


To study the behavior of the internal waves excited by a 
small amplitude disturbance in the presence of a uniform 
oncoming flow, as when a submarine moves through the water at 
constant velocity V, we perform an asymptotic expansion on the 
variables ¢, q, u, v, and y in terms of the small parameter a. 
This is reasonable since the small parameter @ enters the 
bottom boundary condition in an analytic way and since the 
solutions arise by analytic processes, it 1S natural to expect 
the solution to exhibit an analytic dependence on the 
parameter @Q. Since the disturbance is O(@), we shall 
initially suppose that the response is also O(q@). Therefore, 


we introduce the following asymptotic expansions 


etOKG (eee 
GeO yg 2) Clem OP GQiat 1, 
Pieri Uke 0aeee). «FOU aes 37 15) 
We = OLVF OF Vara 


mT AX AY tag Ein Oya tween 


in which we consider the limit process as a > 0. 
Substituting (3.1) into the continuity equation (2.40), we 


obtain 


1S 


0 2 _ 
~z— [ou] i eee + ae pe ee + OlGe) eee = Oe e.2) 


Dividing by a and taking the limit asa-0O, (3.2) simplifies 





EO 
Os a a es (0) (Sram 
where 
D 
W, = oS (3.4) 
DEE 
and 


(355) 





To show that D,/D,t is the O(1) part of the convective 
derivative given in (2.45), we collect only the O(1) terms in 


the expansion for D/Dt given in (3.6) below: 


D 
Dt 


ne a eee (3.6) 
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oe 
ae 
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Zw) 


Substituting (3.1) into the x-component of the momentum 


equations (2.41), we obtain 








p. (2-6) -% (atu,) eG) Cl at OG) ee = 0 (a7) 
Expanding p,(z-¢) in a Taylor Series and using (3.1) yields 
zg) =p (ze pe a toe | Oo) so. (3.8) 
Substituting (3.8) into (3.7), we obtain 
ap,(z) DE me 8 (6 (Wie gael O02) ee meee a0 (3.9) 


Dividing by @a and taking the limit as a > 0, (3.9) becomes 


(2) 2282 BG ao (3.10) 
ae Dit = 
A Similar analysis for the y-component of the momentum 


equations (2.42) yields 


on 





(2) 20Me 0 @ 001) 
ea + = : 
PON el Te ep 
For the z-component of the momentum equations (2.43), we 
substitute (3.1) inte (2243) tevebtarn 
aS il 
ae?" * + S[p,(2=0) =p (2) ) © eq. + OG eo (3.12) 
pt? 8 hs 
Expanding p,(z-f) in a Taylor Series as in (3.8) and 
substituting into (3.12), we obtain 
Zz 
DG Ofis7 
ae’ -2 - 2p, (2) ©, + Gee O(a eer (325085 
pt? 8B x 


Dividing by @ and taking the limit asa->0O, (3.13) simplifies 


Ee 


Po (z) C+ Seon (3.14) 


Notice that the eé* term has been neglected. As previously 


discussed, in order for our theory to work we require & << 1. 
Therefore, we see the se? term is small. The fact that the e’ 


ZZ 


term is small by itself is not justification in dropping this 
term from (3.13). In order to justify the dropping of the e’ 
term, we shall assume that aé’? = o(a@) as a > 0, or that €& = 
o(1) as a> 0. Now as one takes the limit @ > 0 in (3.13), we 
Clearly see that €* > 0. In continuing our derivation of the 
z-component of the momentum equations, we substitute (2.15) in 


for dp,/dz, so that (3.14) becomes 


Ga oe = 0). (aes) 


Next, we consider the surface boundary conditions. We 
begin by substituting (3.1) into the kinematic bottom boundary 


condition (2.46) to obtain 


pee +10 CQ) see at z=af. (3.16) 


(set at z= 0. (34g) 


Note that the O(a) bottom boundary condition is applied at z 
= (0) Substituting (3.1) into the free surface kinematic 


boundary condition (2.47a) results in 


Zo 


GS "Oy + OG") = ee at Za ene (3.18) 


Taking the limit as a > 0, (3.18) becomes 


fo = No at Zl (3.2m 


Again the O(a) free surface condition is applied at z= 1. 
The last boundary condition to consider is the free 
surface dynamic boundary condition (2.47b). Substituting 


(3.1) into (2.47b), we obtain 


p,(1+n) + aBgas O(o7) + -. pasenG (3.20) 


Expanding p,(l+n) in a Taylor Series and using (3.1) gives 


Pp, (1+) = p,(1) + apj(z)n, + O(a?) +... (3.21) 


Substituting (3.21) dnto (3-20) 7 we fine 


P,(1) + op, (4) n, + abies O(el) + yee ene. (3.229 
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Canceling the p,(1) terms, dividing by a, and taking the limit 


as a > 0, (3.22) becomes 


p,(z), + Bq, = 0 at 4 eel (2235) 


but uSing (2.12) simplifies the free surface dynamic boundary 


Genaqition to 


Bq, = 0.1. at Zee yo 


In summary, we see from substituting the asymptotic 
expansion (3.1) into (2.40-2.45), (2.46), and (2.47a,b), that 


the O(a) problem is: 








ee eNO, (oe )) 
Du, 
p,(z) eae pecs Oe (2226) 
DV, 
p,(z) = +g, = 0, (3.27) 
Spt OG = (3.28) 
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where 





with boundary conditions 


(re a at 


6 = No at 


Bq, = PN, ak 


Combining 9(323 2a) sande 25) 


condition at the free surface: 
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gives 


(3 2) 


’ (303 0) 


(3,288 


(3.3224) 


(3.222%) 


the single boundary 


(3 3a 


IV. DEVELOPMENT OF THE TWO-DIMENSIONAL WAVE EQUATION 


A. EIGENVALUE PROBLEM 


In order to solve (3.25-28), (3.31), and (3.33), we look 


memesolutions of the form 


u = Ae! (KX+LyY- ar) gp’ (z) 
Vv = Bei (kX+Ly- ar) g’ (z) 
oO 
(4.1) 
G = Dei (kX + LY - QT) d (z) 
oO 


qd, = Devel (kX + LY - ar) p. gd’ (z) 


where «x = (K,L) 1s the horizontal wave vector, 2 is the 
frequency, and c is the long-wave phase speed relative to the 
basic state at rest. Substituting (4.1) into (3.25), (3.26), 
and (3.27) and solving the system of linear equations, we find 


the dispersion relationship to be 


Queen where wales. (422) 


Furthermore, substituting (4.1) into (3.28), (3.31), and 


(3.33) results in 
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/ / p,N 
(Aro (2) ae —o(z) = 0, (4.3a) 
3 
with boundary conditions 
@(z) = 0 at Z = 108 (4.3b) 
@ (2) = be-O at Ze (4.3c) 


The equations (4.3a-c) formally pose an eigenvalue problem 
where c? is the eigenvalue and ¢ the eigenfunction. Note that 
the dispersion relationship in (4.2) was obtained for wave 
motion in a stratified fluid at rest where V = 0 and D,/D,t = 
d/dt. In the problem we are considering, the dispersion 
relationship would differ from that found in (4.2) but the 
eigenvalue problem would remain the same. Hence, our purpose 
here is to show the existence and character of the solutions 
that define the eigenvalue problem. The solutions have the 
horizontal character of plane waves travelling in the (x,y) 
plane with a vertical structure defined by the eigenfunctions, 
@(z). We shall see that the general solution to (3.25-3.28), 
(3.31), (3.33) can be obtained in terms of a superposition of 
normal modes whose depth dependence is determined by the 


@,.(z), the eigenfunctions of the following eigenvalue problem: 
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N 
Ole | ee rae = 0, (Aga) 
$(z)=0 at Z=0, (4. 4b) 
$,(z) = Bc-d, (Zz) Flee poe il (4.4c) 


Here c,* is the eigenvalue, where c, is the long-wave phase 


speed relative to the basic state at rest. 


B. ORTHOGONALITY 
The eigenfunctions’ ¢,(z) Satisfy the orthogonality 


Geanaition 


oF | Pobs rz dz = oF I,, (4 5) 
0 


where I, is defined by the left side of (4.5) for r = s. To 
prove (4.5), we let A, and A, be eigenvalues with corresponding 
eigenfunctions ¢, and ¢, where A, = c.° and <A, = C,”. The 


differential equations satisfied by the eigenfunctions are 


(4.6a) 


II 
© 


(p,9,, si A, p,N’ p, 


(4.6b) 


I 
© 


(0,6, )' + A, 0,N’ $, 


Zo 


In addition, both ¢, and ¢, satisfy the same set of homogeneous 
boundary conditions defined in (4.4b,c). Multiplying (4.6a) 


by \.¢, and (4.6b) by X,¢, yields 


d,(p,9,. )/o, + A,A,p,N’ >, >, = 0, (4.7a) 
d, (0,6, )'b, + dA, pNP, b, = 0. (4. 7b) 
Subtracting (4.7b) from (4.7a) results in 
A, (9,0, )'b, — A, (0,6, )'¢, = 0. (4.8) 
Through application of the product rule 
h, (050s. ) >, = A-[ (0,6, 6)’ - Doha ba I, (4S 
(4.8) becomes 
(036s OP — Ano, OO Oa pee (4.10) 


Integrating (4.10) yields 
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l 


1 
ae A) [Pods $.4z =. p..$, | - Pods: | (4.11) 


Evaluating the right hand side of (4.11) and remembering #¢,(1) 
= Br,'¢,(1), $,(1) = Bd,'d,(1), and ¢,(0) = ¢,(0) = 0, (4.11) 


becomes 


ee eae po, 6, dz = 0. (4.12) 


ime. =~ A, it follows that 


[j Pots ¢. dz = 0. (4.13) 
0 


In other words, the z derivatives of the eigenfunctions (¢, 
and ¢,) corresponding to different eigenvalues (A, and X,) are 
orthogonal with weight p,. Therefore, the normal modes do 


Satisfy the orthogonality condition in (4.5). 


C. TWO-DIMENSIONAL WAVE EQUATION 

Since the eigenfunctions, ¢,(z), form a complete set, we 
can represent u,, v,, ¢§,, and q, by a generalized Fourier series 
expansion. From (4.1) we know the general form of the 


solutions and hence we put 


sial 


C, = ) A(x, y,t) (2), (4.14a) 


Gs = Po(2)  CzBy Comat (4.14) 
0 
US =) Deca Dee (4.14c) 
0 
Ven =» | Dy a gb) (el 0, ba (4.14d) 
0 


To derive a relationship between A, and B,, we begin with 
C, = Do A,(x, ¥, t) 6, (Zz). (4.15) 
0 


We take the derivative of (4.15) with respect to z, multiply 
by o,c.’¢,(z) and use the orthogonality condition in (4.5) to 


obtain 
al 
I.A, = caf PG. penaae (4.16) 


Integrating (4.16) by parts then gives 


a2 


2 1 
fie As a C5P.6 sz | 7 cof Co (p.,.) dz. (4.17) 
Using (4.4a), (4.17) becomes 


2 | 2 
TAs = CP Ms2| + ip 0.6,N26,dz. (4.18) 


Now, we multiply (4.14b) by ¢, and from the orthogonality 


condition (4.5), we obtain 


1 
I,B, = I, g.o.,dz. (4.19) 
Integrating (4.19) by parts, we obtain 
i 1 
EE os |e Cope (4.20) 
0 0 
Using (3.28), (4.20) becomes 
1 i 
I,B, = a0.) + f o,0,N2o, dz. (4.21) 
Combining (4.18) and (4.21) yields 
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al 1 
CEP. 600s: | + I, B, - qo. | = [ae (A999) 


Substituting (3.33) and (4.4c) for q)(1) and ¢,(1) respeceiaem, 


and remembering that ¢,(0) = £(x,y,t), (4.22) simplifies to 


5+ F, (4.23) 


where 


TF 6 ap Onde OG) alee (4.24) 


Note that F,(x,y,t) 1s proportional to the bottom topography 


Tp Va eae 
Next, we substitute (4.14a,c,d) into the continuity 
equation (3.25) to obtain 








(4.25) 


Multiplying (4.25) by c/o,¢, and applying the orthogonality 


condition (4.5) then gives 
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DA, _ 9D, OE, 


where D,/D,t is given in (3.30). 

To derive the two-dimensional wave equation, we use the x 
and y components of the momentum equations (3.26) and (3.27). 
Bubstituting (4.14b,c) and (4.14b,d) into (3.26) and (3.27) 


respectively, we obtain 








D = oe | ae 
: -5 Dd. | + — S COBO. tea 4.27a 
Po Ds 7 Ps Po x 5 Ds: ( ) 
D = ala 2 
°|-YE + Pow cB. ¢, |= 0: (4.27b) 
P. ID 7 ape Po y ; Or 


Multiplying (4.27a) and (4.27b) by cZ¢, and applying the 


orthogonality condition (4.5) then yields 


POP Cc? HE 0 (4.28a) 
a ae = ; 2 
ie * Ox 
beat ara Bay 0 (4.28b) 
De ‘dy 
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We next take the x derivative of (4.28a) and the y derivative 


of (4.28b) and add the two results to obtain 


_D, fav, _ ae, 
Dt lox ° Oy 





| ee Vioo = — Ue (4.29) 


Substituting (4.23) and (4.26) into (4.29) finally yremee 








D’A 
ae = ¢, V7| A, + Bap (4.30) 
where 
2 2 
ve eee ae (4.31) 
Ox dy’ 
and 
D> F a? 
— - | ae + Ve]. (4.32) 


Equation (4.30) is the inhomogeneous two-dimensional wave 
equation whose solution will be obtained in Chapter V. 
To summarize, we see that substituting (4.14a-d) into 


(3.25-3.27) implies 


a6 








(4.33a) 


(4.33b) 


(43 3c) 


(4233) 


(4.33e) 


(4.22857) 


0 ang Erom (2.7) 


; (4.34a) 


BS tor 
DA, OD, OE 
= + 
DE wo dy 
D,A, 2e7 2 
CNT Ae ieee NG 
Dee. 
where 
Gao) ot Guess tye 
D 0 0 
oO =< ire V : 
Dies Ole Ox 
and 
ee ee US 
Ox? dy’ 
Here, Since f vanishes at t = 0, F,(x,y,0) = 
we have 
vane) re 6G. kay) as eo 
where 


a7 


I,G, = C,p,(0) (0) g(x,y) . (4.34b) 


The initial conditions for (4.14a-d) are that A,, B,, D,, and 


s S 


E, all vanish at t = 0. 
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V. SOLUTION TO THE TWO-DIMENSIONAL WAVE EQUATION 


A. METHOD OF DESCENT 
To solve (4.33c) we first introduce the following change 


of variables: 


C= = Ve, 
ae (5. 1a) 
0 =t 


Derivatives with respect to x, y, t can be expressed as 


derivatives in the new variables ¢, y, o as follows: 


a_ a 
Ox OC’ 

Oy oy’ (Seb) 
cme aC 

Ot 00 0e- 


Substituting (5.1a,b) into (4.33c) yields 


OA, 2, VA , aA, = £ (¢+Ve s (5.2) 
Tre Gee ee a 


a9 


where 


Q? Q? 
£.({+Vo,y,0) = co | ===) ee Ee (Cera (5.3) 
nes y 38° ap is y 


We shall solve (5.2) by beginning with the solution to the 


following three-dimensional wave equation: 


2 2 2 2 

Ou Sey d‘u , du (5.4a) 
Ee Ox? dy’ 02° 
U(X, Vipo 0), = Clee vara, (5.4b) 
u(x, Y¥ 210) = Mie (5.4c) 


Using Kirchoff's formula (see Strauss 1992), the solution to 


(S.4a-c) is 


u(x.,.y,, 2.6 [Youve 2) ds 
ANG ce 


0 aL 
” OE. aor, [Jeter arash 


40 


where 


De | (yz ec eee Veen eect. |. (5.6) 


Hees that the value of u(x,,y,,2,,t,) depends only on the values 
of VW(x,y,z) and $(x,y,z) on the spherical surface D. 
By the method of descent, we now obtain the solution to 


the following two-dimensional wave equation: 





d*u Se eae , Ju, (5.7a) 
at? Ox? dy’ 
(2G, 7 Oy) = Oxy), Cor 0) 
Day 10) Vy) (5a 7 C) 


We regard the solution u(x,y,t) of (5.7a-c) as a solution to 
the three-dimensional problem (5.4a-c) with initial data 


independent of z. By Kirchoff's formula (5.7), we have 





mO.,.01.0). = ee [ ¥onw as 
4nc’t J, 


GSe78)) 


0 aL 
+ ——— | |?(x,y) dS], 
dt resral f 
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where 


+ ya 2? eet |. (5.9) 


This is the solution to (5.7a-c) but it can be simplified. To 
Simplify (5.8), we must express the surface integral as an 


ordinary double integral. It is twice the integral over the 


) 1/2 


top hemisphere z = (c’t’-x’-y’ where (x,y) ranges over the 


region D, = { (x,y) : x? + y’? s c’t? }. By definition 
dS = |R, x R,| dxdy, (5.10a) 


where 


R = eon Co ge) ), Gsaliec. 


and 


(5.10c) 


Substituting (5-10a-e)) intom aece yeas 
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SS ge eee 
HNO, 0 1) xze | | dxdy 


D, Vc? t?-x?-y? 


(Saas) 
daxd 
* ane | J = - 
Ct =x Va 
mois is the solution at the point (0,0,t). Therefore, the 
solution at a general point is 
Me Y57 ©). = ————S=— A dx dy 
oO oO — 
cme -(y-y,)? 
522) 
0 il (x, y) 
+ an dx d 
Ot, a. 4 
D, yc?t.- (x- x,)?-(y-y,)? 
where 
Ia 501d eon ese.) Cen GreN aio csletten |e (Sea) 


Note that (5.12) shows that the value u(x,,y,,t,) depends on 
the values of @#(x,y) and W(x,y) inside the disc defined in 
(5.13). This illustrates the well known fact that Huygens' 
principle fails in two dimensions. We shall use (5.12) in the 
solution to the nonhomogeneous two-dimensional wave equation 


discussed in the next section. 
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B. DUHAMEL'S PRINCIPLE 


We now consider the following nonhomogeneous’ two- 


dimensional wave equation 


07u 9 0° u 07u 
—. —~ Cc ee + —— } = (xX; pes (5.14a) 
dt? dy’ ees 


with homogeneous initial conditions 


Ui cy 0 — Oe (5.14bD) 


Uo V7, 0) = e0' (5. tae) 


By Duhamel's principle, the solution to (5.14a-c) is given by 
ON Sea se) \o U(x, ¥; tasids (5 Ss 

where for each s 2 0, U(x,y,t,s) is the solution of 

542 ae = 0 EORuNe=.sS- (5 Ga 


with initial data prescribed on the plane t = s such that 
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U(x,y,sS,s) = 0, (Sm) 


OA, VS) a9 (Vs). (5. 2oie) 


Since the wave operator is invariant under translations in 


time t (see Strauss 1992), we can shift the time t by a 


constant and still have a solution. We let V(x,y,t,s) = 


U(x,y,t+s,s), so that V(x,y,t,s) is a Solution to 


V. -c*?(V_+V.) =0 for Gee Oe (5.17a) 


with initial data prescribed on the plane t = 0 such that 


V(x, y,0,s) = U(x,y,SsS,s) = 0, (5.615771) 
(> 17 Cc) 


Vale, 0, S) = UX, ) s,s) = Fx, ya ss 


Using Kirchoff's formula (5.12) to obtain a solution to 


(5.17a-c) then yields 


vixy.t.s) = sh { {9 ¥.8)__ ae 
Sy) iekiec= (ej) = (hae 


dy (5-40) 
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where D, = { (x,y) : (x-x)? + (y-y)? s c*t? }. Since our pumpeee 
here is to solve the nonhomogeneous two-dimensional wave 
equation in (5.14a-c) by Duhamel's principle (5.15), we can 
again shift the time t by a constant so that U(x,y,t,s) = 
V(x,y,t-S,Ss). This last expression for U(x,y,t,s) together 


with (5.18) then implies 


U(x,y,t,s) = = OS giles y, 8) dx dy (5.19) 


Vc? (t-s)?-(x-x)?-(y-y)? 


where 


D = [| (x,y) : (x=x)* + (V-y)? = e (e-sie (5.20) 


Now using Duhamel's principle (5.15), we obtain 


SINT dxdy|ds (5.21) 


1 t 
ON Se iE) Paty el? nL Se 
¥ re \J c?(t-s)?-(x-x)?-(y-y)? 


the solution to (5.14a-c) where D, is defined in (5.20). 

In deriving the solution to the nonhomogeneous’ two- 
dimensional wave equation(5.14a-c), we have solved (5.2). We 
change the names of the dependent and independent variables in 


(5.21) to obtain the solution of (5.2). 
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Z (x+Vs,y,s Ss) 
t,y,0) = az, | = dxdy|ds (5.22) 
eto s)?-(x-¢)? 
pieeomme = { (x,y): (x-€)? + (y-y)? = c*(o-s)? }. betting ¢ = 


x - Vt where o = t, (5.22) becomes 


> 

x 

s 
T 


: f a V me ee US Wee) 
| oa sin dx dy |ds 


D> ice (t-s)?- (x- (x-Vt) )?-(y-y)? 


Diesaem = { (x,y) = (x-(x-VE))? + (y-y)? = ci’ (t-s)? }. 

To derive the solution for the nonhomogeneous’ two- 
dimensional wave equation, we began by solving the homogeneous 
three-dimensional wave equation using Kirchoff's formula 
(a5). Next, we used the method of descent to solve the 
homogeneous two-dimensional wave equation by reducing 
Kirchoff's surface integral formula to an ordinary double 
integral (5.12). We developed the solution to the homogeneous 
two-dimensional wave equation in order to solve. the 
nonhomogeneous two-dimensional wave equation by applying 
Duhamel's principle (5.21). This gave the solution of the 


wave equation in (5.2) in the form given above in (5.23). In 


4‘7 


Chapter VI we will examine the solution to (5.23) for long 


time. 
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VI. LONG TIME SOLUTION 


ec 6x VV 

Before looking at the behavior of (5.23) as t > o (where 
the topography f[,(x+Vs,y,s) has a general form), we shall 
first introduce a specific moving point source for the 


topography, so that 


PaaS! a= Ox x ) O4v—-y ) G(s)”. (6.1a) 


Consequently 


f.(x+Vs,y,s) = 6(x+Vs-x,) O(y-y,) 9(s), (6 ae) 


where g is a non-negative function satisfying 


g(0) = 0 and g(s) =1 BOG Se2s. (Oo ane) 


If we can solve (5.23) by substituting (6.1b) in for 
f.(x+Vs,y,s), then we can solve (5.23) for any form of 
topography by applying the principle of superposition. 

Hom solve (5.23), we begin by substituting (6.1b) into 


(5.23) to obtain 
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A ae 


6 (x+Vs-x,) 6(y-y,) g(s) aoe (6.2) 
Fe | Ji - = dx dy |ds. 


We let X = (x+Vs -x,) and Y = (y-y,) and substitute into (6.2) 


tO 7obeain 


A.(x,y,t) = 
2 ———SSSSSSSSSS O (29 0 Y) ats) cingane ae Cae 
271. 
pp en (ese coe) = 
where 
D, =| (X,¥) : [X+V(t-s) -(x-x,)]? + [Y+y,-y]? s Ce (t-s) 0 (6.4) 


In order to solve (6.3) we must determine those s for which 
(6.4) 1s satisfied. Letting 7 = (t-s) for (t-s) = 0, (6may 


Simplifies to 


(c2-V?) 7? + 2V(x-x,)T - [(x-x,)? + (y-y,)’?] 20. (6357 


5) 


Note that X and Y are not included in (6.5) since the 
integration in (6.3) with respect to X and Y makes both X and 


Y become zero. Solving (6.5) for 7 yields 


The values of 7, and 7. will depend on whether c,’? > V’ or c, < 
v’. We shall examine both cases in the following sections. 
me Cc, > Vv 
When the long-wave phase speed is greater than the 
velocity of the submarine, we see from (6.6) that 7, > 0 and 


meeoeO. Factoring (6.5) yields 


el a0) aie) (pais (cay 


where 7, and 7. are defined in (6.6). Since the left hand side 


of (6.7) must be greater than or equal to zero, 7 must be 


greater than or equal to 7, (7 2 7,). Hence, s must be less 
than or equal to (t-7,). Therefore, the upper limit of s in 
(6.3) becomes (t-7,). Now, integrating (6.3) with respect to 


X and Y yields 


oak 


A,(x,y,t) = —__, [""_Ss) & 6.8) 


where 7, and 7. are defined in (6.6). Using (6.1c), we break 


(6.8) into the following parts 


ANA) Y; Gli =2 2 eee (6.9a) 
where 
ee |) SES (6.9b) 
Zc. (Go Va) ((t-7,) -s)((t-7_) -s) 
and 


i, = —— | (6.9C) 


We first look at the behavior of I, as t > ». We see 
from (6.9b) that for long time, I, is O(1/t). Therefore, we 
can neglect this portion of A,(x,y,t). Now in evaluating I, 


for long time, we let 7 = (t-s) so that I, simplifies to 


2) 


=, [| *—S —. (6.10) 


We complete the square of the denominator of (6.10) to obtain 


I, = Se (6.11) 
EMG GeV) oe (7T-a)?-b? 
Here a = (7,4+7.)/2, and b* = (7,-7.)7/4. Letting bu = (7T-a) and 
substituting into (6.11) then yields 
{-s-a 
eet 
ee if aE, (Ga) 
206 (eC, -V) "es, em 
ar ae 


Integrating I,, we find 


Se. 


1/2 


can 26, (c;-V?) 2 ze b* a 
(Gress) 
: We es 
s =e P= 
rey ee - : D ~ + = — | 


For Jong time, 2,915 0 incr). 


In the case where c,’ > V’, 


we see that the long time 
solution for A,(x,y,t) 1s O(ln t). Next we examine the case 
Mopanmigaichey Cla ee 
DT Cee 
When the long-wave phase speed is less than the 
velocity of the submarine, we see from (6.5) that x must be 
greater that x, Since rt 1S greater than or equal to zero. We 


must reexamine (6.6) to see if a positive real value for T 1s 


possible. Rearranging (6.6) for c,” < W yields 


T = ee ee ee (6.14) 
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From (6.14) we see that a real solution for 7 exists only if 


the following inequality holds 


co’ (x-x,)? - (V?-c*) (y-y,)? = 0. (6.15) 


Hence, a solution exists only within the region defined by 


x2x,+ |— -1 vere ey) eT Vee (6.16a) 
Cs 
and 
bie 
Vv? 
x2x,+ |— - 1 (Vie y) bor vo ye (6. 16D) 
Cs 


Since T, and 7. are both positive, we see from (6.7) that 7 
must be between 7, and 7., 1e., (7T.sTsT,). Therefore, s must 
be less than or equal to (t-7T_). 

To solve A,(x,y,t) for c” < V’ as t > », we start with 


(6.3). Integrating (6.3) with respect to X and Y yields 


Sy. 


A,(x,y,t) = : 5 | _ See (6 an 


27C, ( (on (t-7_) -s) 


where 7, and 7. are defined in (6.14). Using (6.1c), we break 


(6.17) into the following parts 


AX V7) =a (6. 2S) 
where 
I, = oa TNSIEes (6.18b) 
2G, (Virem a Were) =Se=7 )) =e) 
and 
ay oe 
Le 3 f SS (6.18c) 
2mC,(C,-V°)"" ¥, (CO=m Ss) (Ge= aa! 
Looking at the behavior of I; as t > », we see from 
(6.18b) “that fits O(a ee Therefore, we can neglect this 
portion of A.(x,y,t). The evaluation of I, is similar to the 


analysis as developed for I,: 
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Wz 


(t-s -a) (t-s-a)? 
= iS : | O o =a | 
ame (V?-e.)\'? b jay 
(6.19) 
1/2 
4 (t_-a) (t_-a)? 
+ ‘aig -1 | 
eme (V2—c:) b b? 
Here a = (7,+7.)/2 and b’? = (17,-7.)?/4. For long time, I, is 
Cin €) . 
In the case where V’* > cZ, we see that the long time 
Hoiieton for A,(x,y,t) 1s O{(ln tt). Next, we examine the 


resonant case where c. = V’. 


B. CC” = V’ 
Previously, we developed the long time solution for the 


* with the interesting result that 


cases when c2 > V? and c2 < V 
in both cases the solution had growth of O(ln t). The reason 
for this growth is the failure of Huygens' principle in two 
dimensions. Though ln(t) tends to infinity for long time, it 
is a function of extremely slow growth, and for the times we 
are interested in, it can be considered as a bounded function. 

In the resonant case, we are concerned with the behavior 


of A,(x,y,t) as t > » for the case when the speed of the 


submarine is equal to the long-wave phase speed (c,’ = V’). 
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To determine if a solution to (6.3) is possible for c? = 


v?, we must reexamine (6.5). From (6.5) we see that a 


solution exists only if s satisfies the following inequality 


Crepe) BOY, - 


EE ——— I ———————— 15OHG V (Xe) > Oe (62209 
2V(%-xX,) 


Integrating (6.3) with respect to X and Y yields 


A(X,Y, G) = 


1 4 g(s) ds (6.21) 
rat a 5 , 
2V(X=x%)) (E=s) — | (Xe) 4 


Here A = [t - ((x-x,)*+(y-y,)?)/(2V(x-x,))]. Using (6.1c), we 


break (6.21) into the following parts 


A Vie) nae (6.22a) 


where 


as 


: Q olerae (6.22b) 
“at 
2V(x-x,) (t-s) -— [ (x-x,)? + (y-y,)7] 
and 
Hi 
r A ae (6. 226) 
a > ; 
8, Y2V(x-X,) (t-s) - [(x-x,)? + (y-y,)7] 
Here A = [t - ((x-x,)7+(ly-y,)7)/ (2V(x-x,))]. 


Looking at the behavior of I; as t > », we see from 
(6.22b) that for long time, I, is O(1/t'”). Therefore we can 
neglect this portion of A,(x,y,t). Integrating I, directly, we 


obtain 


(6a 23) 


From (6.23) we see that I, is O(t'’’) as t > ». Therefore, we 
see that the solution for A,(x,y,t) at resonance grows like 


mie ) for long time. 


5), 


To derive the long time solution (t > ) for A(x)y) wee 
began by introducing a moving point source for the topography, 
f.(x+Vs,y,s), with the knowledge that if a solution was 
attainable by using (6.1c), we could solve (5.23) for any form 
of the topography by the principle of superposition. We first 
looked at the case where c” > VW’ and determined that the long 
time solution for A,(x,y,t) was O(ln t). Next, we looked at 


the case where c,’ < V’ 


and saw that a solution existed only 
within a certain region defined by (6.16a,b). The long time 
solution for A,(x,y,t) in this case was also O(ln t). We 
attribute the O(ln t) growth to the failure of Huygens' 
principle in two dimensions. For the resonant case where c, 
= V’, we found the solution for A,(x,y,t) to be O(t'’) fomemenae 
time. From this analysis, we see in this problem that if a 


Submarine patrols at or near a long-wave phase speed a 


Significant resonant could occur for long times. 
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VII. CONCLUSIONS 


This thesis has presented a theoretical study of the flow 
of a stratified fluid over variable topography as a model of 
the resonant interaction of a submarine's wake with a 
stratified fluid. We considered the case in which a submarine 
patrols at a speed near the long-wave phase speed of an 
internal wave. Internal waves are generated by the flow 
around the submarine's wake. 

The long time solution obtained in our model for the cases 
when c.’ > V’ and c,’? < V’ had growth of O(1n t). (see (6.13) and 
(6.19)). We attribute the O(ln t) growth to the failure of 
Huygens' principle in two dimensions. For the resonant case 
mere c-’ = Vee we found the solution to bese O(t'’?) for long 
time. (see (6.23)). It is important to note that for the 
time scales we are concerned with, both long time solutions 
for the nonresonant cases are bounded and have extremely slow 
growth. Therefore it is the resonant case which is important. 

The growth in time occurring in the resonant case 
indicates that perhaps nonlinear effects, balanced by 
dispersion, could cause significant upstream disturbances. 
From these disturbances, the possibility exists of detecting 


a submarine by the internal waves its wake generates. 
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Further studies in this problem should begin by rescaling 
the asymptotic expansions to balance the leading order 
quadratic nonlinear terms with the topographic forcing. This 
will achieve a balance between nonlinearity and dispersion. 
We expect the solution to this problem to show that 
Significant upstream disturbances could occur when the 
Submarine patrols at a speed at or near a long-wave phase 
speed. 

If internal waves were generated only by underwater moving 
objects then the problem of detecting the submarine by the 
internal waves its wake generates would be made simpler, but 
this is not the case. Internal waves are a natural occurring 
phenomenon in the oceans. This complicates the problem 
because in order to detect the submarine by the internal waves 
that it generates, the ambient internal wave patterns at the 
location of the submarine must be known. The problem of 
discerning ambient internal waves from internal waves 
generated by the submarine's wake is similar to the problem of 
detecting a submarine passively (such as with a tail). The 
sonar operator filters out the ambient noises or background 
noises of the ocean from the incoming signal to detect the 
acoustic signature of the submarine. From further 
understanding of the ocean's naturally occurring internal 
waves, we may one day be able to discern internal waves 


generated by a submarine from those of the ocean. 
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